Nucleus-mediated spin-flip transitions in GaAs quantum dots 
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Spin-flip rates in GaAs quantum dots can be quite slow, thus opening up the possibilities to 
manipulate spin states in the dots. We present here estimations of inelastic spin-flip rates mediated 
by hyperfine interaction with nuclei. Under general assumptions the nucleus mediated rate is pro- 
portional to the phonon relaxation rate for the corresponding non-spin-flip transitions. The rate can 
be accelerated in the vicinity of a singlet-triplet excited states crossing. The small proportionality 
coefficient depends inversely on the number of nuclei in the quantum dot. We compare our results 
with known mechanisms of spin-flip in GaAs quantum dot. 



The electron spin states in bulk semiconductor and 
heterostructures have attracted much attention in recent 
years. Experiments indicate very long spin decoherence 
times and r^igall transition rates between states of dif- 
ferent spincTta. These promising results have motivated 
proposals for information processing based on electron 
spins in quantum dots, ffiifich might lead to a realization 
of a quantum computerou. 

A quantum dot is region where electrons are confined. 
The energy spectrum is discrete, due to, the small size, 
and can display atomic like propertiesotj. Here we will 
consider quantum dots in GaAs-AlGaAs heterostruc- 
tures. The main reasons for studying them are that 
relevant quantum dots are fabricated in such structures 
and GaAs has a peculiar electron and phonon properties 
which are of interest. There are two main types of gate 
controlled dots in these systems, so-called vertical and 
lateral dotstl They are characterized by different trans- 
verse confinement, which is approximately a triangular 
well and a square well for the lateral and vertical dots, 
respectively. 

Manipulation the electron spin in a coherent way re- 
quires that it should be relatively well isolated from 
the surrounding environment. Coupling a quantum dot, 
or any closed quantum system, to its environment can 
cause decoherence and dissipation. One measure of the 
strength of the coupling to the environment are the tran- 
sition rates, or inverse lifetimes, of the quantum dot 
states. Calculations of transition rates between differ- 
ent spin states due to phonon-assisted spin-flip process 
mediated by spin-orbit coupling, which is one possibil- 
ity for spin relaxation, have given surprisingly low rates 
in quantum dotsETllJ. For these calculations it is very 
important that the electron states are discrete, and the 
result differs strongly from that obtained in application 
to 2D extended electron states in GaAs. The same argu- 
ment applies to the phonon scattering mechanisms, since 
certain phonon processes possible in 2D and 3D electron 
systems are not effective in scattering the electron in OD. 
An alternative mechanism of spin relaxation in quantum 
dots is caused by hyperfine coupling of nuclear spins to 
those of electrons. Although the hyperfine interaction 
mediated spin relaxation in donors was considered a long 



time agc0, no analysis has been made, yet, of the hyper- 
fine interaction mediated spin flip processes in quantum 
dots. 

The present publication offers an estimation of the 
scale of hyperfine interaction induced spin relaxation 
rates in GaAs quantum dots and its magnetic field de- 
pendence, the main result presented by the expression 
in Eq. (p^). Since the parameters of hyperfine interac- 
tion between conduction band electrons and underlying 
nuclei in GaAs have been extensively investigatedtJO, 
including the Overhauser effect— and spin-relaxation in 
GaAs/AlGaAs heterostructurealjllfl , we are able now 
to predict the typical time scale for this process in partic- 
ular quantum dot geometries. The rate that we find de- 
pends inversely on the number of nuclei in the dot (which 
can be manipulated by changing the gate voltage) and is 
proportional to the inverse squared exchange splitting in 
the dot (which can be varied by application of an exter- 
nal magnetic field with orientation within the 2D plane 
The following text is organized in two sections: section 
where the transition rates in systems with discrete spec- 
tra are analyzed and section n where the obtained result 
is compared to transition rates provided by the spin-orbit 
coupling mechanism. 



I. MODEL AND ASSUMPTIONS 



The ground state of a quantum dot is a many electron 
singlet \S g ), for sufficiently low magnetic fields. This 
can change at higher magnetic fields. We assume that 
the system is in a regime of magnetic field so that the 
lowest lying states are ordered as shown in Fig. The 
relevant energy scales used in the following analysis are 
given by the energy difference between the triplet state 
(we assume small Zeeman splitting) and the ground state 
c = Et> — Eg, and exchange splitting Sst = Es> — Ej-i 
between the first excited singlet and the triplet, ft is pos- 
sible to inject an electron into an excited state of the dot. 
If this excited state is a triplet, the system may get stuck 
there since a spin-flip is required to cause transitions to 
the ground state. 

The T-point of the conduction band in GaAs is 



f 



mainly composed of s-orbitals, so that the hyperfine in- 
teraction cap, be described by the contact interaction 
HamiltonianE] 



H HF = Aj2Si-IkS(r i -R k ) 



(1) 



where Si (Ik) and (Rk) denote the spin and posi- 
tion the i-th electron (fc-th nuclei). This coupling flips 
the electron spin and simultaneously lowers/raises the z- 
component of a nuclear spin, which mixes spin states and 
provides possibility for the relaxation. 

But the hyperfine interaction alone does not guaran- 
tee that transitions between the above-described states 
occurs, since the nuclear spin flip cannot relax the exces- 
sive initial state energy. (The energy associated with a 
nuclear spin is the nuclear Zeeman, Tiuj n , energy which is 
three orders of magnitude smaller than the electron Zee- 
man energy and the energies related to the orbital degree 
of freedom.) For free electrons, the change in energy ac- 
companying a spin flip caused by the hyperfine scattering 
is compensated by an appropriate change in its kinetic 
energy. In the case of a quantum dot, or any system 
with discrete a energy spectrum, this mechanism is not 
available and no hyperfine induced transitions will occur 
because the energy released by the quantum can not be 
absorbed. Therefore, the spin relaxation process in a dot 
also requires taking into account the electron coupling to 
the lattice vibrations. The excess energy from the quan- 
tum dot can be emitted in the form of a phonon. Since 
the 'bare' electron-phonon interaction, ff p h does not con- 
tain any spin operators and thus does not couple directly 
different spin states, one has to employ second order per- 
turbation theory which results in transitions via virtual 
states. The amplitude of such a transition between the 
triplet state \T') and the ground state \S g ) is 



(T'\S g 
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(T\H ph \t)(t\Hn F \S g ) 
E T , - (E t + hu q ) 

(r'|gHF|s)(g|ff P h|g fl 

Et> - (E a + huj n ) 



(2) 



where fiuj q is the energy of the emitted phonon and 
Tiio n « is the energy changed by raising/lowering a 
nuclear spin. 

It is natural to assume that the exchange splitting is 
smaller than the single-particle level splitting, so that the 
dominating contribution to the amplitude (T'\S g ) comes 
from the term describing the virtual state \s) = \S') , due 
to a small denominator. All other term can be ignored 
and the approximate amplitude takes the form 



(T'\S 9 



(r\H BF \S')(S'\H ph \S g 
Eji/ — Eg/ 



(3) 



The justification for this assumption is that we aim at ob- 
taining estimates of the rates, and including higher states 



would not affect the order of magnitude, even if the ex- 
change splitting is substantial. Note that the phonon and 
nuclear state are not explicitly written in Eq. (^). 

The transition rate from \T') to \S g ) is given by Fermi's 
golden rule 



f HF -ph = ^ £ KT'is,)! 2 ^-^) 



(4) 



where N' q and fj,' are the final phonon and nuclear states 
respectively and E\ and Ef stand for the initial and final 
energies. Inserting Eq. (||) into (^) and averaging over 
initial nuclear states with probability P(fJ-), we obtain 
an approximate equation for the nucleus mediated tran- 
sition rate 

r H F- P h = ^£|(5';^|i/ph|^;^)| 2 5(i?i-i?f) 



E 



p( M )|(r'; / /|g H F|g , ;M)| 2 

(Et> ~E s ,) 2 



(5) 



= rph(£) E {ET ,-E S ,y (6) 

where T p h is the non-spin-flip phonon rate as a function 
of the relaxed energy e = Et< — E$ ■ 

We will approximate the many-body orbital wave- 
functions by symmetric, |\&s), and antisymmetric, \^t), 
Slater determinants corresponding to the singlet and 
triplet states respectively. It is not obvious a priori 
why this approximation is applicable, since the Coulomb 
interaction in few electron quantum dots can be quite 
strongcj. The exact energy levels are very different from 
those obtained by simply adding the single particle en- 
ergies. However the wavefunction will not drastically 
change and especially the matrix elements calculated us- 
ing Slater determinants are comparable to the ones ob- 
tained by using exact ones. The singlet and triplet wave- 
functions can be decomposed into orbital and spin parts: 
\T') = \V T )\T) and \S) = where 

|r) = ^^|i,+i) + ^p|i,-i) + ^|i,o) (7) 

and the hyperfine contribution to the rate becomes 
a iiF = J2P(^\(T'^\H n \S , ; f i')\ 2 



n /J, 

A 2 



Gcorr£(|*l(H fc )| 2 -|*2W| 2 ) 2 (8) 



where \l/i ; 2 are the wavef unctions of the lowest energy 
states and G CO rr contains the nuclear correlation func- 
tions 



(9) 
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where G,, 7 is the symmetric part of the nuclear correla- 
tion tensor, e^ 7K is the totally anti-symmetric tensor and 
the i/jj's are the coefficients in the triplet state expan- 
sion in Eq. (Q). We assume that nuclei are identical and 
non-interacting, which gives G cor r = 1.25 for an isotropic 
nuclear system. 

Let us now introduce the length scales I and zq which 
are the spatial extent of the electron wavefunction in the 
lateral direction and the dot thickness respectively. Let 
C n denote concentration of nuclei with non-zero spin. 
The effective number of nuclei contained within the quan- 
tum dot is 



No 



C n £ 2 zo- 



(10) 



In GaAs N e s ^> 1 and the sum over the nuclei in Eq. (||) 
can be replaced by C„ / d 3 Rk and we define the dimcn- 
sionless quantity 



Tint =FZ0 / tf*Jlfc(|*i(Jifc 



l*2(^)| 5 



(11) 



To relate the hyperfine constant A (which has dimension 
Energy x Volume) to more convenient parameter we note 
that the splitting of spin up and spin down states at max- 
imum nuclear polarization is E n = AC n I, where 7 = 3/2 
is the nuclear spin. Thus, the hyperfine mediated transi- 
tion rate is 



HF-ph 



!>(£) 



E n 
Sst 



Gcorr Tint 
(2/) 2 iVeff 



(12) 



Note that the rate is inversely proportional to the num- 
ber of nuclei N e e in the quantum dot and depends on the 
inverse squace .of Sst, which are both possible to vary in 
experimentsE30. In particular, the singlet-triplet split- 
ting of excited states of a dot can be brought down to 
zero value using magnetic field parallel to the 2D plane 
of the heterostructure, which would accelerate the relax- 
ation process. The nuclear correlation functions in G corr 
may also be naaiainulated by optical orientation of the 
nuclear systemEj'til. 



II. COMPARISON AND ESTIMATES 

We now consider Eq. ([l2]) for a specific quantum dot 
structure. It is assumed that the lateral confinement is 
parabolic and that the total potential can be split into a 
lateral and transverse part. For vertical dots the approx- 
imate transverse wavefunction is 



1/2 



(13) 



where zq is the thickness of the quantum well, i.e. the 
dot thickness. The wavef unctions in the lateral direc- 
tion are the Darwin- Fock solutions <fi n .i(x, y) with radial 



quantum number n and angular momentum I. The sin- 
gle particle states corresponding to (n, I) = (0,0) and 
(0, ±1) are used to construct the Slater determinant for 
a two-electron quantum dot. In the case of these states 
the factor 7i nt in Eq. JTl| ) is then 7j nt = 0.12 for a vertical 
dot and 7i n t = 0.045 for a lateral one. 

One property of the Darwin-Fock solution is the re- 
lation l~ 2 — hujin*/h 2 where w = ^/fi 2 , + /4 is the 
effective confining frequency and u> c = eB/m* the cy- 
clotron frequency. Inserting this and Eq. ( |l0| ) into ( |l2| ) 
the rate for parabolic quantum dots becomes 



THF-ph = T ph (e)hu} - — 
ost 



Eji \ G CO rrTint 



{2If 



h 2 c n z 



(14) 



Spin relaxation due to spin orbit related mechanisms 
in GaAs quantum dots were investigated by Khaetskii 
and Nazarov in Refs. 10 and 0. We will summarize their 



results here for comparison with our hyperfine-phonon 
mechanism. In Ref [l(J it has been found that the dom- 
inating scattering mechanism is due to the absence of 
inversion symmetry. There are three parameters related 
to this mechanism, 



ri = r ph ( £ ) : 



8 /m*/? 2 



r ph (e)6(m*/3 2 )( ff > B B) 



Too 



mo) 4 



r 2 = r ph ( £ ) 



7 {m*(i 2 ){hu) 



24 



E 2 



(15a) 
(15b) 
(15c) 



Here, m*0 2 is determined by the transverse confinement 
and band structure parameters and E z = (p 2 ) / (2m*). 
For vertical dots of thickness zq = 15 nm, then m*/3 2 ~ 
4 x 10~ 3 meV, but one should be cautious when consid- 



ering a different thicknesses, since m* (3 



so the 



rates are sensitive to variations in zq. The transition 



rates in Eqs. (12) and (15a)- (15c) are all proportional to 
the phonon rate T p h evaluated for the same energy dif- 
ference e. It is thus sufficient to compare the only the 
proportionality coefficients. 

Let us now consider for which confining energies the 
different rates are comparable. At zero magnetic field 
Ti = T5 at h£lo w 0.8 meV. The estimated confining en- 
ergies of vertical quantum, dpt s used in experiment are in 
the range (2 - 5.5) meWEiM For those dots T 2 > Ti 
due to the very different dependence on the confinement, 
Ti oc (hflo)~ 3 and T 2 cx hflo. Doing the same for 
Ti and TuF-ph we obtain that those rates are equal at 
Hflo w 4.4 meV. The numerical values used for the hyper.- 
fine rate in Eq..-ffl2|) are the following: E n = 0.13meVt3, 
Sst = 2.3meVEfand a z = 15 nm. Thus, for B = 0T 
and the previously cited experimental values for the con- 
fining energy the dominant transition rate is Eq. (15c). 

For clarity we will give the values of the rates. The 
non-spin-flip rate T p h(£) is given in KN and using the 
values hSlo — 5.5 meV and B = 0Twe get 



3 



r ph w 3.6 x lffs" 1 

THF-ph 

r 2 



2 x 10~ 2 s _1 
1 x 10 2 s -1 . 



(16a) 
(16b) 
(16c) 



The values of level separations are e = 2.7 meV and 
Sst — 2.84 meV are taken from Ref. |2^. 

An application of a magnetic field to the dot may re- 
sult in two effects. The rate T$ becomes larger than r 2 
for magnetic fields around B = IT (hClo = 3meV) to 
B w 5.4 T [KIq = 5.5 meV). More importantly, the ex- 
change splitting Sst can vanish in some cases and the 
hyperfine rate in Eq. (|l2|) will dominate. It is also worth 
noting that in this limit the approximation used in ob- 
taining Eq. (H) becomes very good. Since the rates con- 
sidered here are all linear in the phonon rate the diver- 
gence of T p h at the singlet-triplet transition does not af- 
fect the ratio of the rates. The above estimates were 
focused on vertical dots. To obtain the corresponding re- 
sults for lateral dots the value of 7i n t in Eq. (14) should 
be used. 

In summary, we have calculated the nucleus mediated 
spin-flip transition rate in GaAs quantum dots. The 
comparison of our results to those previously obtained for 
spin-orbit scattering mechanism indicates that the rates 
we obtained here are relatively low, due to the discrete 
spectrum, so that we believe that hyperfine interaction 
would not cause problems for spin-coherent manipulation 
with GaAs quantum dots. Nevertheless, the hyperfine 
rate, which was found to be lower than the spin-orbit 
rates at small magnetic field, may diverge and become 
dominant at certain values of magnetic field correspond- 
ing to the resonance between triplet and singlet excited 
states in the dot. 
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FIG. 1. The lowest lying states of the quantum dot. The 
energy separation of the two singlet states is denoted by e and 
the exchange splitting by Sst- The two rates indicated are 
the phonon rate r p ij and the combined hyperfine and phonon 
rate r H F- P h 
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